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Abstract 

Future wireless networks are expected to be highly heterogeneous with the co-existence of macrocells 
and small cells as well as provide support for device-to-device (D2D) communication. In such muti-tier 
heterogeneous systems centralized radio resource allocation and interference management schemes will 
not be scalable. In this work, we propose an auction-based distributed solution to allocate radio resources 
in a muti-tier heterogeneous network. We provide the bound of achievable data rate and show that the 
complexity of the proposed scheme is linear with number of transmitter nodes and the available resources. 
The signaling issues (e.g., information exchange over control channels) for the proposed distributed 
solution is also discussed. Numerical results show the effectiveness of proposed solution in comparison 
with a centralized resource allocation scheme. 

Index Terms 

Multi-tier cellular networks, device-to-device (D2D) communication, distributed resource allocation, 
auction method. 


I. Introduction 

The future generation (i.e., 5G) of cellular wireless networks are expected to be a mixture of network 
tiers of different sizes, transmit powers, backhaul connections, different radio access technologies that 
are accessed by an unprecedented numbers of heterogeneous wireless devices (1J. The multi-tier cellular 
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architecture of 5G networks, where small cells (e.g., femto and pico cells) arc underlaid on the macrocells 
and also there is provisioning for device-to-device (D2D) communications, is expected to not only increase 
the coverage and capacity of the cell but also improve the broadband user experience within the cell. 
Since uncoordinated power and spectrum allocation for small cells as well as D2D user equipments 
(DUEs) can cause severe interference to the other receiving nodes, efficient resource allocation and 
interference management is one of the fundamental research challenges for such multi-tier heterogeneous 
networks. Due to the nature of the resource allocation problem in multi-tier cellular networks, centralized 
solutions are computationally expensive and also incur huge signaling overhead. Therefore, distributed or 
semi-distributed solutions with low signaling overhead are desirable where the network nodes (such as 
small cell base stations [SBSs] and DUEs) perform resource allocation independently or by the minimal 
assistance of macro base stations (MBSs). 

In this paper, we use the concept of auction from Economics and present a distributed solution to 
the resource allocation problem for a LTE-A based D2D-enabled multi-tier cellular network. The term 
distributed refers to the fact that the SBSs and the DUEs independently determine the allocation with the 
minimal assistance of MBS. The auction approach allows us to develop a polynomial time-complexity 
algorithm which provides near-optimal performance. The auction process evolves with a bidding process, 
in which unassigned agents (e.g., transmitters) raise the cost and bid for the resources simultaneously. 
Once the bids from all the agents are available, the resources are assigned to the highest bidder. 

Despite the fact that there are ongoing research efforts to address the resource allocation problems 
for two-tier networks (e.g., |2j |4|) as well as for D2D communications (e.g., [; 5 ], j6j), the distributed 
solutions for the radio resource allocation problems in a generic D2D-enabled multi-tier scenario has not 
been studied comprehensively in the literature. In |2j, a distributed solution is proposed for a two-tier 
network using evolutionary game. A utility-based distributed resource allocation scheme for small cell 
networks is proposed in (3'j using interference pricing. In |5j and energy and QoS-aware resource 
allocation scheme is proposed for D2D communication. Auction theory has also been used in the context 
of wireless resource allocation problems (e.g., |4], 0 - 0 ). An auction-based subchannel allocation for 
OFDMA and multihop systems is presented in [7j and f8], respectively. However, the power allocation 
is not considered. A resource allocation approach for multi-cell OFDMA system using the concept 
of auction is presented in [9]. The uplink spectrum-sharing resource allocation problem in a sparsely 
deployed femtocell networks is considered in |4] and the authors proposed an auction-based solution. 
However, the works in 0 |9| consider single-tier systems and the works in |2| |4| do not consider D2D 


communication. 
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Fig. 1. A D2D-enabled mutil-tier cellular network. The small cells (such as pico and femto cells) as well as the D2D pairs 
are underlaid within the macrocell by reusing the same radio resources. 


Different from the above works, we apply the auction method to solve the radio resource allocation 
problem in a heterogeneous muti-tier network. We consider a multi-tier network consisting a MBS, a set 
of SBSs (such as pico and femto base stations) and corresponding small cell user equipments (SUEs), 
as well as DUEs. There is a common set of radio resources (e.g., resource blocks [RBs]) available to 
the network tiers (e.g., MBS, SBSs and DUEs). The SUEs and DUEs use the available resources (e.g., 
RB and power level) in an underlay manner as long as the interference caused to the macro tier (e.g., 
macro user equipments [MUEs]) remains below a given threshold. The goal of resource allocation is to 
allocate the available RBs and transmit power levels to the SUEs and DUEs in order to maximize the 
spectral efficiency without causing significant interference to the MUEs. 

The main contribution of this work is a low-complexity decentralized solution to the radio resource 
allocation problem in a multi-tier cellular system. The key feature of the proposed approach is that it incurs 
polynomial time complexity and low overhead for information exchange. We analyze the complexity and 
the optimality of the solution. To this end, we also discuss the applicability of the proposed approach in 
a practical LTE-A based system. 

The rest of this paper is organized as follows. The system model, related assumptions, and the resource 
allocation problem is presented in Section [H] The auction-based distributed solution is developed in 


Section III We present the numerical results in Section IV before we conclude the paper in Section [Vj 
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II. System Model and Problem Formulation 
A. Netw’ork Model and Assumptions 

Let us consider a transmission scenario in a D2D-enabled multi-tier network as shown in Fig. [T] The 
network consists of one MBS and a set of M cellular MUEs, i.e., U m = {1, 2, • • • , M}. There are also D 
D2D pairs and S closed-access SBSs (such as pico and femto base stations) located within the coverage 
area of the MBS. The set of SBSs is denoted by S = {1, 2, • • • S'}. We assume that each SBS serves only 
one SUE during a transmission intervaQ The set of active SUEs is given by W where \U S \ = S. The 
set of D2D pairs is denoted as U A = {1, 2, • • • , i)\. In addition, the d-th element of the sets U At and 
U Ar denotes the transmitter and receiver UE of the D2D pair d 6 U A , respectively. The set of UEs in 
the network is given by U = U'" U W U U A . For ease of presentation, we refer to the SBSs (SUEs) and 
D2D transmitters (receivers) as underlay transmitters (receivers). Therefore, we denote /C T = S U U At 
the set of underlay transmitters (e.g., SBSs and transmitting D2D UEs) and by /C R =U S U U Ar the set 
of underlay receivers (e.g., SUEs and receiving D2D UEs). Hence, K = S + D denotes the total number 
of underlay transmitters/receivers. 

We refer to the small cells and D2D pairs as underlay tier. The SBSs and DUEs are underlaid within 
the macro tier (e.g., MBS and MUEs) since both the macro tier and the underlay tier (e.g., SBSs, SUEs 
and D2D pairs) use the same set M = (1, 2, • • • , A r } of orthogonal RBs. In the considered multi-tier 
network model, each of the network tiers (e.g., macro tier and underlay tier consisting with small cells 
and D2D UEs) has different transmit power, coverage region and specific set of users. 

Each transmitter node in the underlay tier (e.g., SBS and D2D transmitter) selects one RB from the 
available N RBs. In addition, the underlay transmitters are capable of selecting the transmit power from 
a finite set of power levels, i.e., C = (1,2,-•• , L}. Each SBS and D2D transmitter should select a 
suitable RB-power level combination. This RB-power level combination is referred to as transmission 
alignment 

interference caused by the underlay tier to the macro tier. We assume that value of / T[ j is known to the 
underlay transmitters by using the feedback control channels. An underlay transmitter is allowed to use 
the particular transmission alignment as long as the cross-tier interference to the MUEs is within the 
threshold limit. 

We assume that the user association to the base stations (either MBS or SBSs) is completed prior 
to resource allocation. In addition, the potential DUEs are discovered during the D2D session setup by 

'The scheduling of SUEs by the SBSs is not within the scope of this work. 


|. For each RB n G N, there is a predefined threshold l!y U for maximum aggregated 
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transmitting known synchronization or reference signals (i.e., beacons) [10]. According to our system 
model, only one MUE is served on each RB to avoid co-tier interference within the macro tier. However, 
multiple underlay UEs (e.g., SUEs and DUEs) can reuse the same RB to improve the spectrum utilization. 
This reuse causes severe cross-tier interference to the MUEs, and also co-tier interference within the 
underlay tier. Therefore, an efficient resource allocation scheme will be required. 


B. Achievable Data Rate 


(n) 

The MBS transmits to the MUEs using a fixed power p m > 0 for Vn. For each underlay transmitter 


k E /C T , the transmit power over the RBs is determined by the vector Pi = 


’ (1) (U 

Pk A > 


(N) 

>Pk 


i T 


where 


(n) (n) 

p k ' > 0 denotes the transmit power level of transmitter k over RB n. The transmit power p, ’ , Vn must 
be selected from the finite set of power levels C. Note that if the RB n is not allocated to transmitter k, 
the corresponding power variable p k = 0. Since we assume that each underlay transmitter selects only 
one RB, only one element in the power vector Pi is non-zero. 

We denote by g\- the channel gain for the links i and j over RB n. For the SUEs, we denote by up. 
the SUE associated to SBS k E S, and for the DUEs, up- refers to the receiving D2D UE of the D2D 
transmitter k E U At . The received signal-to-interference-plus-noise ratio (SINR) for the any arbitrary 
SUE or D2D receiver, i.e., u k E /C R , k E /C T over RB n is given by 


->) = 
lu k 


( n ) (n) 

9 k, uA 


(n) (n) v-'' V n ) ( n ) i 9 

g m ,UkPm + £ 9 k' uA' + a 

k'eKt T ,k'Ak 


V n ) J . n ) 


(1) 


where <J k l t is the link gain between the SBS and SUE (e.g., up. E ID. k E S ) or the link gain between 
the D2D UEs (e.g., up- E U d . k E U At ), and <imh h is the interference gain between the MBS and the 
UE u k - In the SINR expression <JT|), a 2 = NoWrb where Wrb is the bandwidth corresponding to an 
RB and A)) denotes the thermal noise. Given the SINR, the data rate of the UE up. over RB n can be 
calculated as = W RB log 2 f 1 + liu) • 


C. Formulation of the Resource Allocation Problem 

The objective of radio resource (i.e., RB and transmit power) allocation problem is to obtain the 
assignment of RB and power level (e.g., transmission alignment) for the underlay UEs (e.g., D2D UEs 
and SUEs) that maximizes the achievable sum data rate. The RB and power level allocation indicator for 
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any underlay transmitter k E /C T is denoted by a binary decision variable 


„(«d) 


where 


(n,l) _ 


X', = 


if the transmitter k is trasnmitting over RB n with power level l 
otherwise. 


( 2 ) 


The achievable data rate of an underlay UE Uk with the corresponding transmitter k is written as 

N L 

*« = ££ x^ ,l) W R B log 2 (1 + 7i n) 

n=11=1 

The aggregated interference experienced on RB n is given by 

K L 

(n,l) in) (ra) 




( 3 ) 


( 4 ) 


k= 1 1=1 


where m* k = argrnax g^ n , Vm E W". The concept of reference user [11] is adopted to calculate the 

m ' 

aggregated interference using (|4j). For any RB n, the interference caused by the underlay transmitter 


k is determined by the highest gains between the transmitter k and MUEs, e.g., the MUE m* k who is 
the mostly affected UE by the transmitter k. Satisfying the interference constraints considering the link 
gain corresponding to the reference user will also satisfy the interference constraints for other MUEs. 
As mentioned earlier, an underlay transmitter is allowed to use a particular transmission alignment only 
when it does not violate the interference threshold to the MUEs, i.e., E"- 1 < , Vn. 


Let X = 


x 


(id) 




)Xi j ■ ■ ■ > x^ 


(■ N,L ) 


> X K 


(N,L) 


i T 


denote the resource (e.g., transmission align¬ 


ment) allocation vector. Mathematically, the resource allocation problem can be expressed as follows: 


(PI) 

K N L 

4 n ’° ^RB log 2 (l + 7^ } ) 

fc=l n= 1 1=1 



subject to X]2 x fc l ’ 0 ^?m;Pfc n) < 4 h> Vn 

k=11=1 

(5a) 


N L 

££4”’‘> = i. v* 

n=l1=1 

(5b) 


x^ E {0, 1}, Vk,n,l 

(5c) 


where 


ry(n) — 

* Ufa 


( n) (n) 

9 k,u k Pk 


(■ n) ( n ) . V' vd ( n,l ') (n) (n) n 

9m.u k Pm + E E x k- 9k' u k Pk' + 


( 6 ) 


k'&K T ,l '=1 
k'^k 


The objective of the resource allocation problem PQ] is to maximize the data rate of the SUEs and 


DUEs. With the constraint in (5a), the aggregated interference caused to the MUEs by the underlay 
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transmitters on each RB is limited by a predefined threshold. The constraint in (5b I indicates that the 
number of RB selected by each underlay transmitter should be one and each transmitter can only select one 
power level at each RB. The binary indicator variable for transmission alignment selection is represented 
by the constraint in (Del. Note that the decision variable = 1 implies that p f ” 1 = l. 


Corollary 1. The resource allocation problem P[l] is a combinatorial non-convex non-linear optimization 
problem. The complexity to solve the above problem using exhaustive search is of O (^(N L) K j and the 
centralized solution is strongly NP-hard for the large values of K, N, L. 


Due to mathematical intractability of solving the above resource allocation problem, in the following 
we present a distributed solution using tools from the auction theory. The distributed solution is developed 
under the assumption that the system is feasible, i.e., given the resources and parameters (e.g., size of the 
network, interference thresholds etc.), it is possible to obtain an allocation that satisfies all the constraints 
of the original optimization problem. 


III. Distributed Solution Using Auction Method 

The resource allocation using auction is based on the bidding procedure, where the agents (i.e., underlay 
transmitters) bid for the resources (e.g., RB and power level). The transmitters select the bid for the 
resources based on the costs (e.g., the interference caused to the MUEs) of using the resource. The desired 
assignment relies on the appropriate selection of the bids. The unassigned transmitters simultaneously 
raise the cost of using resource and bid for the resources. When the bids from all the transmitters are 
available, the resources are assigned to the highest bidder. 

In an auction-based assignment model, every resource j associated with a cost Cj and each agent i 
can obtain benefit Bij from the resource j. The net value (e.g., utility) that an agent i can obtain from 
resource j is given by B r] — Cj. For an equilibrium assignment, every agent i should be assigned with 
resource j' such that the condition B l j> — Cj' > max {/i, ? — Cj } — e is satisfied for all the agents, where 


e > 0 indicates the parameter related to the minimum bid requirement [121. In the following we utilize 
the concept of auction in order to obtain the distributed solution of the resource allocation problem. 


A. Utility and Cost Function 


i(n,Z) a_ (n) 


Let us introduce the parameter I,,,/ ' = 7 u fc j p (n) = j that denotes the achievable SINR of the UE u/. over 

(n) (n) !i~ 

RB n using power level l (e.g., p k ; = I) where p u f is given by (16 1). We express the data rate as a 






function of SINR. In particular, let 31 = Wrb log 2 ^1 + denotes the achievable data rate 

for transmitter k over RB n using power level l. 

The utility of an underlay transmitter for a particular transmission alignment is determined by two 
factors, i.e., the achievable data rate for a given RB and power level combination, and an additional cost 
function that represents the aggregated interference caused to the MUEs on that RB. In particular, the 
utility of the underlay transmitter k for a given RB n and power level l is given by 


4 n ’° 




(7) 


where the operator [-] + = max{0, •} and v \, are the biasing factors which can be selected based on 
which network tier (i.e., macro or underlay tier) should be given priority for resource allocation (2]. Note 

represents the cost (e.g., interference caused by underlay transmitters 
to the MUE) of using the RB n. In particular, when the transmitter k is transmitting with power level l, 
the cost of using RB n can be represented by 


that the term ^ (/< 


n) _ r( n ) 


TH 


in ( 



= *7 


&t l + 


fc'erV'=i 

k'^k 


( n,l ') (n) (n) 

! fc' 91^1 



(8) 


Let the parameter cj?’^ = 


,( n -0 


and accordingly the cost = 


cost term we can represent <71 as where 

data rate achieved by transmitter k using resource {n. /}. 


0 only if I <n) < Using the 
= uiJl is the weighted 


B. Update of Cost and Bidder Information 

Let denotes the local bidding information available to transmitter h for the resource { a, / }. In 

the beginning of the auction procedure, at any time slot t, each underlay transmitter updates the cost as 

- !)} • W 

In addition, the information of highest bidder of the resource {n,l} is obtained by 

b^ l \t) = b^: l \t-i) do) 

where k* = argrnax | (t — 1), cj^’ l \t — 1)|. When the cost of {n, 1} is greater than previous 

k'&K T ,k'^k *• J 

iteration and the transmitter k is not the highest bidder, the transmitter needs to select a new transmission 
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Algorithm 1 Auction method for any underlay transmitter k 


Input: Parameters from previous iteration: an assignment X(t — 1) = [xi(t— 1)] , aggregated interference 

I(t — 1) = \f n \t — 1)1 , cost of using resources C(t — 1) = \ci n ’ l \t — 1)1 and the highest bidders of the 

L Jvn L Jvt.n.i 

resources 23(t - 1) = [93 fc (t - l)]y fc where <8 fc (-) = [bi n ’ !) (-)l 

L J Vn,Z 

Output: The allocation variable Xfe (t) = , updated costs Cfc(t) = C^ n,!) (f) , and bidding information *Bfc(t) = 

L J Wn,l L J Vn,Z 

^ at curren ); iteration t for the transmitter k. 

L J Vn,Z 

1: Initialize Xfc(t) := 0. 

2: For all the resources {n, 1} £ JV x £ update the cost and highest bidder using (|9]l and i 10 1 , respectively. 

3: {n, 1} := Non-zero entry in Xk(t — 1). /* Resource assigned in previous iteration */ 


4: if c£' l) (t) > - 1) and f k then 




(n,l) , 


5: 

6 : 

7: 

8 : 

9: 

10 : 

11 : 

12 : 

13: 


{h,l} := argmax iijj." ,l ft). 

{n' ,1' } GAT X C 

T 1 "-* := g^J + I (n) . /* Estimate interference level */ 

if < I™ then 


Set x 


(ft,0 


:= 1. 


Update the highest bidder for the resource {n, 2} as bf’ l \t) := fc. 


Increase the cost Cjf'^ (t) = — 1) + Ak(t — 1) where Afc(t — 1) is given by 0- 


else 


Xfe(f) := Xk(f — 1). /* Keep the assignment unchanged */ 

end if 
14: else 

15: Xfc(t) := Xfc(t — 1). /* Keep the assignment unchanged */ 

16: end if 


alignment, say, {h,l}. The transmitter also increases the cost of the new resource {h. / } as C^" J ' (t) = 
— 1) + Afc(t — 1), where A k (t — 1) is given by 


A k (t-1) = max U, 


(n',V) 


(Cli¬ 


max 


it 1 




{n f xC 


(f - 1) + e. 


( 11 ) 


r«*'V + _-n_ max if («'.*')/ 

{n',Z'}eA/'x£ 

the maximum and the second to the maximum utility value. In the case when the transmitter k does 


The term max llT ’(t — 1) — max 111“ ,l; (t — 1) physically represents the difference between 

{n',Z'}eA/’x£ ~ k r.. 


not prefer to be assigned with a new resource, the allocation from the previous iteration will remain 

T 


unchanged, i.e., x k (t) = x k (t — 1), where x/ ; = 


A n d) 


. Wn,l 


C. Algorithm for Resource Allocation 

We outline the auction-based resource allocation approach in Algorithm [2] Each transmitter locally 
executes Algorithm |T| and obtains a temporary allocation. After execution of Algorithm [l| each underlay 
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Algorithm 2 Auction-based resource allocation 
Phase I: Initialization 

1: Estimate the channel state information from the previous time slot. 

2: Each underlay transmitter k £ fC T randomly selects a transmission alignment and reports to the MBS. 

3: MBS broadcasts the assignment of all transmitters, aggregated interference of each RB, the costs and the highest bidders 
using pilot signals. 

4: Initialize number of iterations t := 1. 

Phase II: Update 

5: while X(t) f X(f — 1) and f < T max do 

6: Each underlay transmitter k £ JC T locally runs the Algorithm[l] and reports the assignment Xfc(f), the cost C k(t) and 

the bidding information 23fc(t) to the MBS. 

7: MBS calculates the aggregated interference vector 1(f); forms the vectors X(t), 23(f), and C(f); and broadcast to the 

underlay transmitters. 

8: Update f t + 1. 

9: end while 

Phase III: Transmission 

10: Use the resources (e.g., the RB and power levels) allocated in the final stage of update phase for data transmission. 


transmitter k reports to the MBS the local information, e.g., choices for the resources, x/,. = 


„fiM) 


i T 


. Vn,Z 


Once the information (e.g., output parameters from Algorithm [T]) from all the transmitters are available 
to the MBS, the necessary parameters (e.g., input arguments required by Algorithm [I]) are calculated 
and broadcast by the MBS. Algorithm [l] repeats in an iterative manner until the allocation variable X = 


I x fc]vfc = 




, 


„(!.£) 


, 


„{N,L) 




1 T 


for two successive iterations remains unchanged. 


D. Optimality, Complexity, and Signaling over Control Channels 

Proposition 1. The sum data rate obtained by the distributed auction algorithm is within Ke of the 
optimal solution. 


Proof: Refer to Appendix [A] 


Since the proposed scheme satisfies the allocation constraints (|5b|)-(|5cl), and also maintains the inter¬ 
ference threshold given by ( [5a] ), the solution obtained by the auction algorithm gives a lower bound of 
original resource allocation problem P[T| As shown in the following proposition, when the allocation 
remain unchanged for at least T > 2 consecutive iterations, the complexity of the proposed solution is 
linear with number of underlay transmitters and the available resources. 
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Proposition 2. The auction algorithm converges to a fixed allocation with the number of iterations of 


O (: TKNLT) where T = 


inuxB;" ' 1 min B 

k,n,l k,n,l 


(n,l) 


Proof: Refer to Appendix [B| ■ 

In the following, we discuss the applicability of the proposed method in a practical system. In the 


proposed solution, the MBS only needs the allocation vector X and the channel gain vector G = 

1 T j 

to calculate the aggregated interference vector, e.g., I = ” ] Vn - Also note that at each 


(n) 

9 k, ml 


k J Vfc,n 

iteration, the MBS and the transmitters need to exchange the allocation vector X, information about 
highest bidders 23. the cost vector C to the underlay transmitters. However, in LTE-A based systems, 
these information can easily be incorporated into the standard control channel messages. For example, 
the MBS can broadcast these information vectors using physical downlink control channel (PDCCH). 
Likewise, each of the underlay transmitters k can inform the MBS the local assignment preferences x/,., 
updated costs C+ and bidding information 23/- using physical hybrid-ARQ indicator channel (PHICH). 


IV. Numerical Results 

A. Simulation Setup 

We develop a MATLAB-based simulator and observe the performance of our proposed approach using 
simulations. We simulate a 300 m x 300 m area where the MBS is located in the center of the cell 
and M = 6 MUEs are uniformly distributed within the cell radius. The SBSs and SUEs are uniformly 
distributed within the macro cell and small cell radius, respectively. The DUEs are located according to 


the clustered distribution model 1131. We obtain the minimum bid increment parameter e by trial and 


error and set it to 100. We choose number of RBs N = 6, biasing factors = U 2 = U and assume the 
interference threshold to be —70 dBm for all the RBs. 

For modeling the propagation channel, we consider distance-dependent path-loss and shadow fading; 
furthermore, the channel is assumed to experience Rayleigh fading. In particular, we use realistic 3GPP 


propagation environment presented in [14]. For example, the following path-loss equation is used to 
estimate path-loss between SBSs and SUEs as well as to the MUEs and DUEs: 'ips(^)\dB\ = 38.46 + 
201og(T) + £ ss + 101og(?). The path-loss between the MBS and MUEs as well as DUEs and SUEs is 
estimated as V'M {V)[dB] = 15.3+ 40 log(k) + £ sm + 101og(Y) + £ w . Similarly, the direct link gain between 
DUEs is given by ipr)(£)[dB] = 148 + 401og(0.001f) + Gd + 101og(c) + £ w . In the path-loss equations, 
£ is the distance between nodes in meter, +. s , £ sm , account for shadow fading and are modeled as 
a log-normal random variables, g is an exponentially distributed random variable which represents the 
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TABLE I 

Simulation Parameters 


Parameter 

Values 

Cell layout 

Isolated, hexagonal grid 

Transmission bandwidth 

1.08 MHz 

Macro and small cell radius 

300 m and 30 m 

Distance between D2D peers 

15 meter 

Number of SUE per SBS 

1 

Transmit power of MBS 

43 dBm 

Standard deviation for shadowing: 


Macrocell and D2D pairs 

8 dB 

Small cells 

4 dB 

Outdoor wall loss 

30 dB 

Noise power spectral density 

-174 dBm/Hz 


Rayleigh fading channel power gain, and £ w denotes outdoor wall loss. Numerical results are averaged 
over different independent spatial network realizations and channel gains. The key simulation parameters 
and assumptions are listed in Table [T] 


B. Results 

The convergence behavior of the proposed solution is depicted in Fig. [2j We can find that the data 
rate of the network becomes stable within 100 iterations. In order to observe the convergence behavior 
in different network density, we vary the number of underlay nodes (e.g., SBSs, SUEs and DUEs) 
and plot the empiri cal cu mulative distribution function (CDF) of the number of iterations required for 

The empirical CDF is defined as F r (j) = 1 where r is the total 


convergence in Fig 


2(b) 


i= 1 


number of simulation observations, Q is the number of iterations required for convergence at the v-th 


simulation observation, and j represents the x-axis values in Fig. 2(b) The indicator function Ij.j outputs 
1 if the condition [•] is satisfied and 0 otherwise. As mentioned in Proposition [2j the convergence to 
the fixed allocation depends on number of underlay transmitters. When the number of network nodes 
increases, the number of iterations required for convergence increases. This is because, the underlay 
transmitters need to execute the algorithm more times in order to obtain the updated bidding and cost 
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information, and hence to find the fixed allocation. However, even in a moderately dense network with 30 
underlay nodes (e.g., S + D = 9 + 6 = 15), the solution converges to a fixed rate within 100 iterations. 



(a) 



(b) 

Fig. 2. (a) The convergence of proposed solution (for L = 3, e.g., C = {3,5,7} dBm, S = 6, D = 4), and (b) Empirical 

CDF of the number of iterations required for convergence using similar power levels as those for Fig. |2(a)| 

To evaluate the performance of the proposed algorithm, in Fig. [3] we compare the the data rate with 
the upper bound obtained from the solution of the original optimization problem P[l] We obtain the 
upper bound (e.g., maximum achievable data rate) using exhaustive search. We examine the sum rate 
of the network in a period of 50 time slots and average the results over different spatial realizations. 
We measure the sum rate efficiency as r? = B prop where R „ rop and f? max denote the sum data rate 
obtained from proposed scheme and maximum achievable sum rate, respectively. From this figure and 
as mentioned in Proposition [l] it can be observed that the data rate of the proposed solution is within 
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Fig. 3. Achievable sum rate of the proposed scheme and the upper bound. We consider S = 3, D = 2, and L = {3, 5} dBm. 


Ke of the maximum achievable data rate. Recall that the original resource allocation problem is an 
NP-hard problem and the computational complexity of exhaustive search to achieve the upper bound 
is of O ((NL) j. The proposed scheme provides data rates which are close to 80% of the maximum 
achievable data rate, however, with significantly less computational complexity and signaling overhead. 

V. Conclusion 

We have presented a polynomial time-complexity distributed solution approach for the future hetero¬ 
geneous multi-tier cellular mobile communication systems. We analyze the optimality and complexity 
the solution, and discuss about the applicability of the proposed scheme in a practical system. Numerical 
results have shown that the proposed solution performs close the upper bound of achievable data rate with 
significantly less complexity and minimal involvement of the central controller node. As an extension of 
this work, capturing the dynamics of misbehaving and selfish nodes in future dense networks using the 
other auction-based method (e.g., truthful auction) could be an interesting area of research. 


Appendix 


A. Proof of Proposition [7] 


We construct the proof by using an approach similar to that presented in |15|. The data rate obtained 
by any arbitrary assignment using auction method will satisfy the following condition: 


E e< "'°+E { „,5^ 

k= 1 {n,l}£j\f xC k =1 


(A.l) 
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where C tn J ^ = maxC^ and R Uk is given by (31. 


k'£lC T 

The inequality in (|A. 1 [) is satisfied since 


K N L 


E = EEE<A 


(n,l) n (n,l) 


{n,l}£Af xC 


and 


K 


k=ln= 11=1 

K N L 


a EEE 4 ”'VT' ) - 

fc=i L 


fc=ln=l1=1 




K N L 


Let the variable A = max ^ v\R Uk = ^ W'rtB log 2 (1 + 7«”j denotes the maximum 

X k =1 k=ln=ll=l ' ' 

achievable weighted sum rate. In addition, let the variable D* be defined as 

D* 4 min { y C^ l) + J2 
{npeJVxC {n,l}£MxC k =1 

Hence by definition we can write A* < D*. Since the final assignment and the set of costs are at 


max 

(n,l}sA/’x£ 


{ 5 (n,l) _^( n ,l ) }j ( A . 2 ) 


equilibrium, for any underlay transmitter k, the condition 

N L 

V V 4 n ’° (B ( k n,l) - C (n ’°) > max {s[ n ’° - 

k v k ' - {n,QeMxC l k J 

n— 1/=1 

will hold. Consequently, we can obtain the following inequality: 

K N L 

d * < y f y y + max {y ,o - 

k {npeAfxC 1 k 

k= 1 n= 11=1 


— e 


K N L 


<EEE4 n ’ 1, BT i, +^ 


lc=l n=l 1=1 
K 


< y + iTe < A* + Xe. 


(A.3) 


fe=i 


Since A* < D*, the data rate achieved by the auction algorithm is within Ke of the optimal data rate 
A* and the proof follows. 


B. Proof of Proposition [2] 

From mu it can be shown that, in the worst case the total number of iterations in which a resource 


(e.g., transmission alignment) can receive a bid is no more than T = 


1" ’ 1J — min P ^" ’ 1 ' 

k.n.l k.n.l 


. Since each 


bid requires O ( NL ) iterations, and each iteration involves a bid by a single transmitter, the total number 
of iterations in Algorithm [2] is of O (KNLY). For convergence, the allocation variable X needs to be 
unchanged for at least T >2 consecutive iterations. Hence, the overall running time of the algorithm is 
of O (TKNLT ). It is worth noting that for any transmitter node k G /C r , the complexity of the auction 
process given by Algorithm [l] is linear with the number of resources for each of the iterations. 
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